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ABSTRACT 

The galaxy power spectrum contains information on the growth of structure, the 
growth rate through redshift space distortions, and the cosmic expansion through 
baryon acoustic oscillation features. We study the ability of two proposed experiments, 
BigBOSS and JDEM-PS, to test the cosmological model and general relativity. We 
quantify the latter result in terms of the gravitational growth index 7, whose value 
in general relativity is 7 ~ 0.55. Significant deviations from this value could indicate 
new physics beyond the standard model of cosmology. The results show that Big- 
BOSS (JDEM-PS) would be capable of measuring 7 with an uncertainty 17(7) = 0.043 
(0.054), which tightens to (7(7) = 0.031 (0.038) if we include Stage III data priors, 
marginalizing over neutrino mass, time varying dark energy equation of state, and 
other parameters. For all dark energy parameters and related figures of merit the two 
experiments give comparable results. We also carry out some studies of the influence 
of redshift range, resolution, treatment of nonlinearitics, and bias evolution to enable 
further improvement. 
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1 INTRODUCTION 

Surveys of large-scale structure in the universe provide a 
rich resource for testing our understanding of cosmology. 
Future surveys will cover nearly the full sky to redshifts 
far deeper than are currently studied, mapping out some 
10 billion years of history. The great statistical power and 
leverage from depth will allow detailed examination of the 
cosmological framework by carrying out a simultaneous fit of 
a substantial suite of relevant parameters. One particularly 
attractive prospect is the capability to put to the test the 
predictions of Einstein gravity for the growth of structure 
and its consistency with the cosmic expansion history. 

We consider next-generation surveys mapping the dis- 
tribution of galaxies in three dimensions to redshifts of order 

2 = 2. A goal of this study is to determine the capabilities 
of such surveys. In particular we aim to estimate realistic 
constraints from a global parameter fit on the gravitational 
growth index 7, which can characterize deviations from gen- 
eral relativity. The second goal is to examine how the survey 
characteristics such as redshift range, resolution, and galaxy 
selection affect those capabilities. 

In Sec. [2] we review the formalism for extracting cosmo- 
logical information from galaxy correlation measurements 
in terms of the matter power spectrum, and discuss the 
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anisotropic distortion due to measuring in redshift space 
(rather than position space). We discuss the relevant set 
of cosmological parameters in Sec. [3] and their influence on 
the matter power spectrum. The results are analyzed with 
emphasis on the role of degeneracies between factors that 
influence growth, including the gravitational growth index, 
the dark energy equation of state, and neutrino mass. In 
Sec. [4] we turn to astrophysical and survey characteristics 
and analyze the effect of the bias level of the selected galaxy 
populations, the form of the small-scale velocity damping, 
the spectroscopic survey redshift resolution, and the redshift 
range of the survey. This allows quantitative comparison of 
the capabilities of next-generation (Stage IV) experiments 
from both ground and space, as well as nearer term (Stage 
III) experiments. We conclude in Sec. [5] with a summary 
of the prospects for testing the standard cosmology and re- 
vealing clues to dark energy or the breakdown of Einstein 
gravity. 



2 METHODOLOGY 

The future dark energy experiments considered in this pa- 
per aim at measuring galaxy positions in three dimensions 
to study baryon acoustic oscillations and other aspects of 
the matter power spectrum including its evolution through 
the growth of structure. The matter power spectrum con- 
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tains important cosmological information through its evolv- 
ing amplitude, its shape including the turnover reflecting 
the transition from radiation to matter domination and the 
suppression due to massive neutrino free streaming, and the 
baryon acoustic oscillation features serving as a standard 
ruler. 

One aspect of particular interest is the distorted, 
anisotropic mapping between the real space density field and 
the mea surements in redshift space, c aused by peculiar ve- 
locities (|Kaiserl Il987l : Irlamiltonl [l998h . This redshift space 
distortion has attracted recent attention as a possible tech- 
nique for dete c ting deviation s from genera l relat ivity (see 
IPeeblesl (|2002l ); iLinderl (|200Sl ); iGuzzo et al.l (|2008l ) for early 
work) as it depends on the relation between the density and 
velocity fields, which can be altered by modifying the grav- 
itational theory. 

Thus the observed galaxy power spectrum contains sev- 
eral types of cosmological information. The autocorrelation 
function £(r) is defined as the excess probability of finding 
masses at a separation r: 



dV = p(l + S m )dV, 
dVu = p 2 (l+^(r))dV 1 dV 2 



(1) 

(2) 



where p is the mean mass density and 8 m = (p — p)/p is 
the density contrast. The mass power spectrum is then the 
Fourier transform of the autocorrelation function: 



P(k) 



d 3 rC(r)e* r , 



(3) 



with k the wave vector. Due to spatial isotropy, only the 
magnitude k will enter. 

We do not observe the power spectrum in real space, 
however, but obtain the radial position through redshift 
measurements, convolving the real distance with additional 
redshifts due to peculiar velocities. This leads to the redshift- 
space power spectrum P gainin g an angular dependence 
through the linear Kaiser factor (|Kaiserl [19871 ) multiplying 
the isotropic, real space mass power spectrum P(k): 

P(k, t J.) = (b + f t /) 2 P(k), (4) 

where [i is the cosine of the angle that k makes with the line 
of sight. For notational simplicity, we suppress the tilde from 
now on. We work in the linear regime, where the continuity 
equation between the galaxy peculiar velocity fiel d and the 

alaxy mass overdensity is linear (see for example iHamiltonl 

199S 



The dimensionless growth rate / is given by 



/ = 



d\nD 
din a ' 



(5) 



where a is the scale factor, and D(a) is the growth factor, 
i.e. the amplitude <5 m (k, a) oc D(a) or P(k) oc D 2 (a). We 
also need to take into account that galaxies, not directly 
mass density, are observed. The bias b relates the galaxy 
overdensity S g to the total mass overdensity through S g — 
b5 m . 

By looking at the angular dependence of the power spec- 
trum at each k, 



P(k,p) oc a s (b + ) 



alb 2 + 2a%bfn 2 + of/ 2 : /i , (6) 



2 2 2 2 2 

in principle fit for b a 8 , bfa 8 and / a 8 , hence allowing us to 
measure b and / provided we have an appropriate measure- 
ment of a 8 . This is challenging in practice due to noise. An- 
other possible route to separating out th e bias involves the 
use o f higher order correlation functions l|Scoccimarro et all 
GUI). 

Although we have three measurable quantities (the 
three coefficients of the fourth order polynomial in Eq. [6j 
and three unknowns, we cannot determine all of them be- 
cause the second is the geometric mean of the other two. 
This is because we work in the linear regime and general rel- 
ativity, where the galaxy density and peculiar velocity fields 
are perfectly correlated. But should one of these hypothe- 
sis be relaxed (as in modified gravity models or with non- 
linearities e.g. Finger-of-God effects) , we need to introduce 
the correlation c oefficient between the fields (|White et al.l 
l2009l ; IUzanll2009l ) 



r(k) = 



Pg V {k) 



(7) 



where the subscript g denotes the galaxy density field, 
and v the divergence of the peculiar velocity field. Ideally 
this correlation would be predicted by the physical theory 
ijDesiacques et ai~1l2009t ): allowing r instead to b e completely 
free s ignificantly degrades the constraints on / (|White et al.l 
I2009D . We do not consider this situation further in this arti- 
cle, instead assuming the standard correlation of unity, since 
we restrict our analysis to the linear regime and many classes 
of gravity theory maintain the correlation in this regime. 

To incorporate a measure of the sensitivity to the grav- 
it y theory we us e the gravitational growth index formalism 
which parameterizes the growth factor as 



lt y tneory we us 
of lLinderl l|2005h 



da' 



D(a) = a exp / [fi m (a' ) 7 - 1] — , 



/ = n m { a y 



where 



£3 m (a) 



(8) 



(9) 



(10) 



where as is the normalization of the power spectrum, we can 



is the ratio of matter density to the total energy density at 
scale size a = (1 + z) _1 . The summation runs over all the 
different components of the universe: matter, dark energy, 
curvature and radiation. The gravitational growth index 7 
will be a parameter of key interest. It can distinguis h other 
theori e s from Einst e in gr a vity - (see for exa mple iLinderl 
|2005l ); lLinder et all (|2007h : lGuzzo et al.l (|2008h ')- The merit 
of a large scale structure survey in terms of its gravitational 
probative power may be conveniently quantified by the un- 
certainty a(y). The F igure of Merit Science Working Group 
jAlbrecht et al.ll2009T ) found that for the suite of future Stage 
III experiments, expected to be completed before the pro- 
posed Joint Dark Energy Mission (JDEM) program, the an- 
ticipated uncertainty is (7(7) = 0.21. 

The standard technique for making suc h parameter es- 
timat ion predictions is the Fisher matrix (|Tegmark et al.l 
I1997D . For a survey covering a volume Vo where the mean 
galaxy number density is n, the element of the Fisher matrix 
for parameters pi and Pj is obtained as an integral over the 
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space of modes k (lTegmarkl ll997h. by: 



p. . — 

~ 2(2tt)3 



d s k ( nP{k,») 



d\nPd\nP 



1 + nP(k, (j,) J dpi dp 



(11) 



The accessible modes are weight ed due to shot noise 1/n 
according to an effective volume (jFeldman et alj|l994l ): 



V„(k,n) = V a 



nP(k, fi) 



(12) 



1 + nP(k, /Lt) j 

The constraint leverage comes mostly from regions where 
nP{k,n) > 1, that is V e « V . 

In order to avoid the uncertainties associated with treat- 
ment of non-linearities, we truncate the Fisher matrix in- 
tegral at a maximum value k+. We take k+ — 0.1 fr/Mpc, 
which is the scale where departures from line ar theory begin 
to bec ome significant (see, e.g., the analysis of lPercival et alj 
l|2009h ). See Sec. l4.2l for a further investigation of non-linear 
effects. 

The information about 7 comes from two different parts 
of the power spectrum. The real space, isotropic part, cor- 
responding to no redshift space distortions, or = (obser- 
vations transverse to the line of sight) in the linear regime, 
is proportional to the growth factor squared: 

\2 



P±{k) = b z P{k) cc D(aY 



(13) 



Note that surveys lacking sufficient redshift resolution are 
only sensitive to the transver se modes due to smear ing along 
the line of sight (see, e.g., iPadmanabhanl (|200&f 0. Using 
Eq. ((5}, the information carried by this part involves 

^^ = 2 rn m (a') 7 lnO m (a')^. (14) 
"7 Jo a 

The redshift space distortions in the power spectrum give 
further information through the parameter /, which with 
Eqs. (0 and flSJ) reads 

f = n m ( a y. (15) 

Therefore, if we define the anisotropic part alone as 

P aniso (fc, At) = 2&/V + /V, (16) 
it carries information on 7 through 



dlnP aniso . . f2 m (a) 7 lnf2 m (a)p 2 

— m 1 Lm \€b) + " 



26 + n m (a)V 



(17) 



This factor gives a sense of the information from the redshift 
distortions. 

Because the measurements become noisier when subdi- 
vided into angular bins, and because a substantial majority 
of the info rmation resides in the spherically averaged power 
spectrum rtShoii et al.|[2009T ) , analyses frequently use the one 
dimensional, spherically averaged power spectrum 



P sph (fc) =P{k) ( b 2 + hf + lf 



(18) 



This incorporates information from both the original 
isotropic power spectrum and the redshift distortion 
anisotropies, and may be most familiar in te rms of the Dy oc 
\D 2 A /H (z)] 1/3 oc (fcifc||r 1/3 factor of lEisenstein et"ai1 
(20051 ) ■ In particular, the sensitivity to 7 arises from 



<91nP E 



spli 



dlnD 2 

97 



+ 



[lQbQ. m {a)-< +6» m (a) 27 ]lnn m (a 

27 



15fo 2 + 106fi m (a)T + 30n(a) 



BigBOSS 


LRG a 


EL 




— 1 


1 — 2 




Z4UUU 




n (n/Mpc) 


3.4 X 10~ 4 


3.4 X 10 -4 


b 


1.7 


0.8 - 1.2 


R 


^ 2300 


> 2300 


JDEM-PS 


LRG a 


EL 


2 range 


- 0.7 


0.7 - 2 


f^skv (deg 2 ) 


10000 


20000 


n (h/Mpc) 3 


3.4 x 10~ 4 


19.5 x 10" 4 


b 


1.7 


0.8 - 1.2 


R 


ks 2000 


^ 200 



(19) 



Table 1. Survey specifications for the Stage IV experiments Big- 
BOSS and JDEM-PS. "Uses northern hemispher e (10000 deg 2 ) 
LRG z = - 0.7 from BOSS jSchlegel et al.ll2009h . 



In Sec. O we will investigate the relative importance of the 
transverse, anisotropic, spherically averaged, as well as full 
versions of the power spectrum for constraints on the grav- 
itational growth index and other parameters. 



3 PARAMETER CONSTRAINTS 

The constraints on 7 expected from nearer term (Stage III) 
surveys are not that informative, as mentio ned, with (7(7) = 
0.21 compared to a difference A 7 = 0.13 ()Lue et all 120041 : 
iLinderlfeoOSl ; iLinder et al. 2007) between general relativity 
and DGP gravity (|Dvali et all I2OO0I : iDeffavet et"ai] 120021 ) 
for example. We therefore turn to Stage IV experiments and 
assess their potential for a more accurate test of the standard 
cosmological model. 

We consider two ve rsions of Stage IV p ower spectrum 
experiments: BigBOSS (|Schlegel et al.ll2009h is a proposed 
ground-based wide fi eld spectroscopic survey and JDEM-PS 
(|Gehrels et al.l |2009| ) is a proposed space-based wide field 
grism survey. Both aim at measuring the three dimensional 
spatial distribution of galaxies to study baryon acoustic os- 
cillations and the growth of structure. Bot h experiments 
woul d use the Stage III experiment BOSS (|Schlegel et al.l 
120091 ). detecting luminous red galaxies (LRG) out to z = 0.7, 
as a springboard to higher redshifts. BigBOSS would ex- 
tend mapping of LRG out to z = 1 and to the southern sky 
and both experiments would supplement LRG with different 

classes of emis sion line galaxi e s (EL ) out to z ~ 2. 

Fo llowin g ISchlegel et all (|2009h ; iGehrels et all (|2009h : 

ISlosarl (|2009l ). we give in Table [1] the redshift range, survey 
solid angle Q s k y , expected target galaxy bias factors bLRG 
and bEL, mean galaxy number density n, and wavelength 
resolution R — A/ A A of the spectrographs to be used (so 
the redshift resolution a z = 5z/(l + 2:) = R^ 1 ). We consider 
variations in redshift, number density, and bias in Sec. [3] 

To calculate the power spectrum as a function of red- 
shift and cosmological par ameters we u sed the Boltzmann 
equation code CMBeasy (|Doranl l2005h . Using two sided 
derivatives together with convergence tests we can accu- 
rately calculate the sensitivity derivatives with respect to 
each parameter. These then enter into the Fisher matrix cal- 
culations of the parameter estimation, taking into account 
the correlations between parameters. The data points are 
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taken to be the power spectrum evaluated at the centers of 
10 (or 11) redshift bins from z = — 2, i.e. at Zi = 0.2i + 0.1. 
For JDEM-PS, we divide the bin containing z — 0.7 into 
two pieces: z = [0.6,0.7] using LRG and z = [0.7,0.8] using 
EL. 

The parameter set involves 9 parameters. Note that 
when testing the gravitational framework, i.e. exploring 
beyond-Einstein gravity through quantitative estimation of 
7, it is crucial to include all parameters that could act in 
a similar manner on the growth and growth rate. There- 
fore we include a time varying dark energy equation of state 
w(a) = wo+w a (l — a) and massive neutrinos. The parameter 
list, and the fiducial value around which the Fisher matrix 
expands, is 

(i) 7 = 0.55, gravitational growth index 

(ii) bhRG, the bias for LRG (see Table [1} 
(hi) &el, the bias for EL (see Table [T} 

(iv) Qde = 0.744, dark energy density today 

(v) Q.v — 0.002, massive neutrino energy density today 

(vi) Wj, = Q, a h 2 = 0.0227, reduced baryon energy density today 

(vii) h = flo/(100 km/s/Mpc) = 0.719, reduced Hubble con- 
stant 

(viii) wq = —0.99, dark energy equation of state today 
(ix) w a — 0, dark energy equation of state time variation 

The val ues for Qde, uib and h are WMAP-5 best fit 
parameters (|Hinshaw et alj [2009). Note that the fiducial 
7 = 0.55 is the value predicted by General Relativity for 
ACDM (and is quite insensitive to the dark energy equation 
of state); the fiducial wq = —0.99 is taken to avoid issues 
of stepping over w = — 1. Dark energy perturbations are 
included in CMBeasy. We assume there is no spatial curva- 
ture. In the remainder of this section we take the fiducial 
bEL = 0.8, and we will investigate the effect of a different 
fiducial in the next section. Note that the neutrino energy 
density fraction is related to the sum of the neutrino masses 
by fi„/i 2 = ^^7n„/94 eV. For a reasonable current upper 

bound < 0.3 eV jSeliak et al.l l2009h . this implies 

fi„ < 0.006. We take 0„ = 0.002, or^m„= 0.1 eV as the 
fiducial. 

Adding together the information from the redshift slices 
independently (note this is not generally a good approxima- 
tion for slices thinner than our Az = 0.2), we obtain the 
full Fisher matrix. We do not explicitly add any CMB infor- 
mation (except later when adding Stage III Fisher matrices, 
which assume Planck data). 

Concentrating on testing the gravitational growth in- 
dex, we now explore in more detail what affects the con- 
straints on 7 using information only from the galaxy power 
spectrum. The constraints are computed to be 



o"(7)sigBOSS = 0.043 
<t(7)jdem-ps = 0.054 . 



(20) 
(21) 



The importance of including dark energy properties, neu- 
trino masses, and other cosmological parameters in the pa- 
rameter estimation is highlighted by the much tighter con- 
straints obtained if we neglect their influence, including 
only 7 itself and the galaxy biases. In this case we ob- 
tain overly optimistic estimates: o"(7)Bi g BOSS = 0.0096 and 
c(7)jdem-ps = 0.0078. Thus, taking into account the corre- 



^BigBOSS 



/0.9954\ 
0.9943 
0.9911 
0.9933 
0.9993 
0.9893 
0.9990 
0.9997 

Vo.9996/ 



^JDEM-PS 



/0.9970\ 
0.9608 
0.9960 
0.9908 
0.9994 
0.9895 
0.9988 
0.9997 

Vo.9996/ 



Table 4. Vectors of the global correlation coefficients for the 
parameters (7, bLRG, bEL, Qde, &u, <^>bi h, wo, w a ) for BigBOSS 
and JDEM-PS. 



lations with other cosmological parameters is essential. The 
correlation matrices for the two experiments are shown in 
Tables [2] and [3] we have replaced the unit diagonal with the 
uncertainties en on each parameter. 

To obtain an overall view of how tightly correlated a 
parameter is with the other variables, we employ the global 
correlation coefficient - the largest correlation of that pa- 
rameter with any linear combination of all other parameters. 
This is given by 



n = Jl- 



Fu {F- 



(22) 



We show those vectors in Table [4] Note the high degree 
of correlation, indicating the importance of crosschecks by 
other data and techniques. 

Examining the marginalized parameter estimations 
along the diagonals of Tables [2] and [31 we see that as expected 
the power spectrum information is especially strong in con- 
straining e and h. One can determine at the ~ 10% level 
the growth index 7 and present equation of state wo, while 
w a and Q. v have uncertainties of order unity. The growth 
index and equation of state parameters estimation is similar 
for the two experiments: 7 = 0.55±0.043, wo = — 0.99±0.16 
and w a = ± 0.47 for the ground-based BigBOSS and 
7 = 0.55±0.054, w = -0.99±0.14 and w a = 0±0.37 for the 
space-based JDEM. We find the usual high anti-correlation 
between wo and w a , and a strong correlation between 7 and 
(wo,w a ). 

Regarding the neutrino mass parameter, neutrino os- 
cillation experiments indicate tha t neutrinos do have mass 
(|Maltoni et al.ll200Sl ; iKavserllioOsT ) , but this is not always in- 
cluded in parameter estimation despite its correlations. We 
demonstrate the effect of neglecting this ingredient, finding 
that it gives overly optimistic constraints on 7 by a factor of 
three to four. The results in Table illustrate the influence 
of neutrinos in three ways, including their mass as a free 
parameter, including their mass but fixing its value, and ne- 
glecting their mass so they act as a relativistic species. At the 
level of neutrino energy density used as fiducial, Q. v — 0.002, 
and over the range k < 0.1 /i/Mpc used for the power spec- 
trum, the parameter value does not strongly affect determi- 
nation of 7 and is mainly degenerate with the bias parame- 
ters. However it is crucial to include neutrino mass because 
the difference between treating them as relativistic vs. non- 
relativistic energy density is still important. 

It is interesting to explore where the main information 
on the gravitational growth index comes from between the 
transverse and anisotropic parts (Eqs. [131 and [T6[), and to 
compare with the spherically averaged case (Eg. I18p . Note 
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/0.043 


-0.50 


-0.40 


0.45 


-0.30 


0.05 


0.31 


-0.93 


0.88 \ 


-0.50 


0.021 


0.96 


-0.70 


0.72 


-0.09 


-0.34 


0.55 


-0.71 


-0.40 


0.96 


0.0099 


-0.70 


0.69 


-0.08 


-0.35 


0.41 


-0.58 


0.45 


-0.70 


-0.70 


0.0039 


-0.12 


-0.10 


0.84 


-0.61 


0.71 


-0.30 


0.72 


0.69 


-0.12 


0.0021 


-0.29 


0.33 


0.18 


-0.33 


0.05 


-0.09 


-0.08 


-0.10 


-0.29 


0.00049 


0.01 


0.00 


0.01 


0.31 


-0.34 


-0.35 


0.84 


0.33 


0.01 


0.0092 


-0.51 


0.53 


-0.93 


0.55 


0.41 


-0.61 


0.18 


0.00 


-0.51 


0.16 


-0.97 


V 0.88 


-0.71 


-0.58 


0.71 


-0.33 


0.01 


0.53 


-0.97 


0.47 / 



Table 2. BigBOSS correlation matrix for the parameters (7, feiijc feei,, HdSi f^, wi,, h, wo, w a ). The off-diagonal elements are r;j = 
Cij/y/ CaCjj while the diagonal elements have been replaced with <r; = \J Cu in bold. 



/0.054 


-0.11 


-0.33 


0.33 


-0.22 


0.01 


0.19 


-0.93 


0.82 \ 


-0.11 


0.018 


0.91 


-0.57 


0.75 


-0.05 


-0.10 


0.26 


-0.52 


-0.33 


0.91 


0.0080 


-0.62 


0.80 


-0.05 


-0.11 


0.41 


-0.63 


0.33 


-0.57 


-0.62 


0.0028 


-0.14 


-0.14 


0.76 


-0.57 


0.70 


-0.22 


0.75 


0.80 


-0.14 


0.0019 


-0.23 


0.43 


0.17 


-0.37 


0.01 


-0.05 


-0.05 


-0.14 


-0.23 


0.00039 


0.02 


0.03 


-0.02 


0.19 


-0.10 


-0.11 


0.76 


0.43 


0.02 


0.0066 


-0.44 


0.44 


-0.93 


0.26 


0.41 


-0.57 


0.17 


0.03 


-0.44 


0.14 


-0.95 


V 0.82 


-0.52 


-0.63 


0.70 


-0.37 


-0.02 


0.44 


-0.95 


0.37/ 



Table 3. JDEM-PS correlation matrix for the parameters (7, b^nQ, bEL, Ode, ^u, h, wo, w a ). The off-diagonal elements are r;j = 
Cij I y/CuCjj while the diagonal elements have been replaced with Oi = \J Cu in bold. 



Case BigBOSS JDEM-PS 

Massive neutrinos, Q„ free 0.043 0.054 
Massive neutrinos, Q v fixed 0.042 0.053 
Relativistic neutrinos 0.014 0.013 

Table 5. Gravitational growth index uncertainty 17(7) under dif- 
ferent treatments of neutrino mass. 

Case BigBOSS JDEM-PS 

Transverse ( - Ea.ll3l 0.126 0.128 

Spherically averaged fEq.[T8t 0.081 0.065 
Full (Eq.gJ 0.043 0.054 

Table 6. Gravitational growth index uncertainty (7(7) using dif- 
ferent parts of the power spectrum. 

that we have defined the anisotropic part to isolate the red- 
shift distortion, imagining one could remove all shape (k) 
dependence and only focus on the angular dependence. This 
seems unrealistic and is only included as a toy model to high- 
light the 7 influence on the growth rate /; the constraints on 
7 become 0.023 for BigBOSS and 0.021 for JDEM-PS (note 
that the parameter space is much reduced, with the baryon 
density, neutrino density, and h not entering). Table [6] shows 
the more realistic parts. 

Note the full power spectrum with redshift space distor- 
tions has the greatest information on the growth index, with 
a factor 2 better constraints than the spherically averaged 
power spectrum and a factor 3 better than the transverse 
(zero redshift distortion or 2D) modes, for the BigBOSS 
case. BigBOSS achieves these improvements due in large 
part to its high resolution that lets it probe the redshift 
distortions more successfully. 

Finally, a Stage IV power spectrum experiment will 
not exist in isolation. Previous experiments, using several 
methods, will be carried out and the complementarity be- 
tween methods offers leverage to tighten the cosmology con- 



BigBOSS +m JDEM-PS+in 

(7(7) 0.031 0.038 

<r(w ) 0.105 0.094 

a{w a ) 0.340 0.289 

Table 7. Gravitational growth index and dark energy equation of 
state uncertainties provided by each of the Stage IV experiments 
in conjunction with Stage III. 



straints. To study the impact of Stage III priors on the pa- 
rameters we_usetJie_Stage_III matrix given by the FoMSWG 
website (|Albrecht et al.|[2009h (without double counting the 
BOSS information), rotated into the (wo,w a ) basis. Sum- 
ming the Fisher matrices of our analysis and of Stage III, 
we extract the constraints on cosmology shown in Table [7] 
The complementarity of the other methods (supernova 
distances, CMB power spectra, and weak lensing shear) from 
Stage III in breaking degeneracies tightens the constraints 
on 7 produced by BigBOSS and JDEM-PS by a factor of 
1.4. Stage IV experiments using these techniques will fur- 
ther reduce the uncertainties on 7, either directly or indi- 
rectly through constraining other, correlated cosmological 
parameters. 

BigBOSS from the ground and JDEM-PS from space 
appear comparable in their cosmology reach. For the 
marginalized uncertainties, BigBOSS does better on the 
gravitational growth index 7 by a factor 1.26 while JDEM- 
PS does better on the equation of state time variation w a 
by a factor 1.27. We exhibit the joint 68% confidence con- 
tours in Fig. [T] where we see that the JDEM-PS contours 
are slightly fatter, having an overall area 1.23 times the Big- 
BOSS constraints. Treating the inverse area of the parame- 
ter estimation contours as a figure of merit (FoM), Table [S] 
lists the ratios of FoM's for the BigBOSS plus Stage III and 
JDEM-PS plus Stage III experiments. 
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Figure 1. lcr joint confidence contours for the gravitational 
growth index 7 and equation of state time variation w a , marginal- 
izing over the other parameters, are plotted for BigBOSS and 
JDEM-PS with and without Stage III information. 



BigBOSS/JDEM-PS BigBOSS+m/JDEM-PS+m 



7,n DE 0.93 0.99 

■y,w 1.16 1.20 

7,w a 1.21 1.23 

w ,w a 0.88 0.86 



Table 8. The ratios of the figures of merit (inverse areas) are 
given for various parameter spaces listed in the first column. The 
second column shows the ratios for the Stage IV experiments 
alone; the third column includes Stage III information for each of 
them. 



4 SURVEY CHARACTERISTICS 

In this section we investigate the influence of different survey 
parameters in the determination of the gravitational growth 
index 7. We discuss the influence of the redshift resolution 
cr z = R~ , the model for non-linear redshift distortions (i.e. 
the small scale velocities appearing in the Finger-of-God ef- 
fect), the uncertainty in the bias parameter bEL shown in 
Table [T] and in particular the survey redshift range and de- 
sign. We will use the information coming from the full power 
spectrum as defined in Eq. (Q. To clarify the effects we do 
not include information from Stage III experiments. 



4.1 Redshift Resolution 

The effect of the uncertainty in the redshift measurement is 
incorporated by including a Gaussian suppression factor in 
the power spectrum in fc-space: 

Pda mp (fc, n) = P(k, fj.) e - fe V^ 2 /«M 2 _ (23) 



R = a- 1 fc* j0 (h/Mpc) BigBOSS JDEM-PS 

20 0.0067 0.110 0.124 

200 0.067 0.044 0.054 

00 00 0.043 0.054 

Table 9. Impact of the resolution R on (7(7). 

We include this factor in the Poissonian noise factor enter- 
ing the effective volume, Eq. (|12fl . but do not vary it with 
cosmology. 

A simple rule of thumb can be derived for the mini- 
mal resolution to achieve in order to neglect the influence 
of redshift uncertainties. Given that we truncate the inte- 
gral defined in the Fisher matrix at fc+ = 0.1/i/Mpc to 
exclude non-linear redshift distortions, this resolution effect 
will start to be significant when fc* = H(z)/ (ca z ) w k+ . This 
yields a z w 0.003 or R « 300. This is near the JDEM-PS 
minimal resolution, but this estimate is for the worst-case 
scenario (k = k+, [J. = 1, H(z) = Hq), so redshift uncer- 
tainties should not be an issue for JDEM-PS/BigBOSS. For 
experiments with larger redshift measurement uncertainties, 
however, the effect on cosmology determination can be sig- 
nificant as shown for the full numerical computations in Ta- 
ble n 

As expected, the BigBOSS/JDEM-PS values fc* = 
0.067 — 0.67/i/Mpc are sufficient for the cosmology estima- 
tion in our studies. However, two issues must also be kept in 
mind: including information from k > 0.1 /i/Mpc would in- 
crease the resolution requirements, and high resolution plays 
a key role in cleanly selecting the galaxy populations, e.g. 
avoiding line confusion in emission line galaxies. 

4.2 Non-linearities 

We can examine how a better understanding of the transi- 
tion to the non-linear part of the power spectrum could lead 
to an improvement in determining a(-y). Instead of trun- 
cating the integral in Eq. flllj l entering the Fisher matrix 
at k+, we can choose to implement a streaming mo del rep- 
resent ing Fingers of God effects (see for example IPeeblesI 
(1980)). where we integrate over all k but multiply the power 
spectrum by a damping factor, either Lorentzian or Gaus- 
sian. This is supposed to model an exponential or Gaussian 
probability distribution function for the peculiar velocities 
of galaxies. We investigate the three forms of the small-scale 
velocity damping factors: 

Cutoff: P n i{k,n) = P{k, n) @(k+ - k) (24) 
Gaussian: P n i(k,n) = P(k, y) e ' {k/h+)2 ^ (25) 

Lorentzian: P nl (k,y) = - -%^2 2 ( 26 ) 

1 + (k/k + ) 2 fi' ! 

where Q is the Heaviside function. We use k+ = 0.1 /i/Mpc 
for all cases. Table [lOl shows the effects on the determination 
of the gravitational growth index. 

We see that the statistical uncertainty on 7 is largest if 
we simply cut out all translinear information, by about a fac- 
tor 2. Thus we have adopted the most conservative method 
to predict 17(7); the information might not be completely 
lost on translinear scales, but only attenuated by Finger- 
of-God effects. Adopting a Gaussian or Lorentzian damping 
model allows extraction of some information, with the choice 
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Case 



BigBOSS JDEM-PS 



Cutoff 
Gaussian 
Lorentzian 



0.043 
0.024 
0.019 



0.054 
0.026 
0.021 



Table 10. The impact of different models for the translincar 
damping due to peculiar velocities on the gravitational growth 
index estimation 17(7). To restrict to the translinear scales we 
further truncate the power spectrum integral at k = lh/Mpc. 



Populations 



BigBOSS JDEM-PS 



LRG, b LRG = 1.7 0.067 0.115 

EL, b EE = 0.8 0.574 0.187 

EL, b EL = 1.2 0.503 0.197 

LRG + EL, b EL = 0.8 0.043 0.054 

LRG + EL, b EL = 1.2 0.042 0.053 

Table 11. Impact of the redshift range and the associated two 
different populations on (7(7). The top three lines consider a single 
population and its redshift range from Table [T] while the bottom 
two lines combine both populations and their redshift ranges. The 
second vs. third, and fourth vs. fifth, lines examine the effect of 
different values for b E L- 



of model affecting the results at the ~ 25% level. However, 
an exponential or Gaussian probability distribution function 
for the streaming model is still not completely accurate, and 
along with the reduced statistical uncertainty on 7 could 
come a systematic bias. Thus we retain the conservative, 
cut off method. T aking into account a halo model, for exam- 
ple [Tinkei] |2003), could allow a more detailed investigation 
of the proper treatment of the translinear regime. Issues of 
nonlinear bias could also arise beyond the k+ = 0.1 ft/Mpc 
adopted in this paper. 



4.3 Redshift Range and Survey Design 

While the survey volume due to the solid angle f2 s k y sim- 
ply scales the parameter estimation as 17(7) oc 1/ y/ fi s k y in 
the statistical treatment without priors, the influence of red- 
shift range is more complex and interesting. Since the galaxy 
population used also depends on redshift we simultaneously 
investigate the influence of the galaxy bias values. 

Table [11] shows the results for considering the popula- 
tions, and their associated redshift ranges, one at a time and 
also in combination with different values (0.8 vs. 1.2) for the 
emission line gal axy population bias. 

As found in iLinderl (|2008l ) , most of the constraint on 
7 comes from the redshift range 2 < 1, which mostly cor- 
responds to the LRG population. The reason is simple: the 
cosmological information on 7 enters the power spectrum 
through the factor r2 m (z) 7 , so at higher redshifts where 
Q m (z) is closer to 1, the sensitivity to 7 decreases. The value 
of the EL bias adopted does not have a significant effect, 
especially when in combination with the low redshift, LRG 
sample. Furthermore, note that the EL only case for JDEM- 
PS, which includes all the information from JDEM-PS itself 
and none of the data to be provided by BOSS, only deter- 
mines (7(7) « 0.2, even though the sample extends down to 
z = 0.7. For JDEM-PS, the BOSS data enable an improve- 
ment of almost a factor 4 in the growth index parameter 
determination. 



These consequences of redshift range raise an important 
question: what is the science reach of the BigBOSS survey 
if the EL sample is shifted from 2 = 1 — 2 to 2 = 0. 7 — 1.7? 
This not only changes the redshift range of the EL sam- 
ple information but creates an overlap between LRG and 
EL information. The g eneralization of Eq. (Ill p to multi- 
ple g alaxy populations (|McDonald et all 120081 ; White et alj 
2009) reads: 



dpi 



dpj 



(27) 



where X and Y are indices describing pairs of galaxy pop- 
ulations, and Cxy is the covariance matrix of the power 
spectra. Adapting the BigBOSS specifications from Table [T] 
by shifting the EL sample to 2 = 0.7— 1.7 retains the science 
leverage and in fact delivers a mild improvement of 8%: 



BigBOSS standard 
BigBOSS z EL = 0.7-1.7 



(7(7) = 0.043 
(7(7) = 0.040 



(28) 
(29) 



Moreover, a redshift maximum of 1.7 reduces the tech- 
nical complexity of the data acquisition and analysis, greatly 
ameliorating issues of line confusion and reduced signal-to- 
noise that occur over 2 = 1.7 — 2. (Note that for 2 > 2 Ly-a 
enters the spectral range and the issues again disappear.) 
The overlap of LRG and EL populations with very different 
biases in the same redshift range 2 = 0.7 — 1.0 also offers the 
possi bility of crosscorrelation and reduction of sample vari- 
ance (|McDonald et all 120081 ) . Thus, these results motivate 
shifting the EL redshift range to 2 = 0.7 — 1.7, achieving 
(7(7) = 0.040 (and 0.030 with Stage III information). 



4.4 Galaxy Samples 

The values of the galaxy number densities and biases listed 
in Table[T]come from the references given. While it is beyond 
the scope of this paper to do detailed survey design, we 
can explore whether some variations in the adopted values 
matter. 

As we have seen in the previous subsection, a change in 
the constant bias of the ELG population from 0.8 to 1.2 has 
a 2% effect on determining 7. We now consider an evolv ing 
model for bias. Motivated by IPadmanabhan et all (|2006h we 
ta ke bLRG = fei + 0-42 wit h fiducial b\ = 1.6, and motivated 
bv lSumivoshi et all <|2009h we take b E L = 62 + (z - 0.7)/2.6 
with fiducial b% = \. 

The constraints on the dark energy parameters 7, u>o, 
w a improve by 4%, 2%, 6% for BigBOSS and degrade by 
5%, 11%, 10% for JDEM-PS. These changes are due to 
altered covariances between the bias parameters and the 
dark energy parameters, involving an interplay between the 
nP factor in the effective volume and the Fisher sensitivity 
9 In P/dbi. Note that the latter quantity goes as 2/(6i + ///), 
so an increased bias decreases the Fisher element. However, 
increasing the bias increases the effective volume through 
raising nP. In the BigBOSS case, this second factor is more 
than sufficient to compensate for the reduced sensitivity. 
However, JDEM-PS has such a high galaxy number den- 
sity that the change in nP has little effect on the effective 
volume, leaving only the reduced sensitivity. Updating Ta- 
ble [8] for the evolving bias case, Table [121 shows further gains 
in the figures of merit for BigBOSS relative to JDEM-PS. 
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BigBOSS/JDEM-PS BigBOSS+m/JDEM-PS+ni 



or 0.031 and 0.038 respectively when combined with nearer- 
term, Stage III experiments. This represents nearly an order 
of magnitude improvement over Stage III knowledge. 

We have also studied the survey characteristics and con- 
firm that the power spectrum at redshifts z < 1 has strong 
leverage. This makes the luminous red galaxy component of 
the survey quite important. Furthermore, our results demon- 
strate that shifting the redshift range of the emission line 
galaxy survey of BigBOSS from z = 1 — 2 to z = 0.7— 1.7 can 
improve the constraints, while adding benefits such as re- 
duced technical complexity and line confusion and increased 
signal-to-noise and the ability to crosscorrelate galaxy pop- 
ulations of different biases. 

Lyman-a forest spectra from BigBOSS quasars at z > 
2, which we have neglected, will further advance the deter- 
mination of cosmological parameters. 

The prospects for testing standard cosmology and in 
particular general relativity are promising. Improved un- 
derstanding of the translinear density regime and velocities 
would further extend the number of usable power spectrum 
modes, while complementarity with other Stage IV exper- 
iments utilizing supernova distances, CMB measurements, 
and weak lensing data would give powerful leverage on both 
the gravitational growth index and other cosmological pa- 
rameters. The capability of probing beyond-Einstein gravity 
opens up a new window for our understanding of cosmic ac- 
celeration and fundamental physics. 
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7,H DE 1.19 1.12 

7, wq 1.45 1.41 

7,tu Q 1.52 1.42 

w ,w a 1.13 1.07 

Table 12. As Table [8] but using an evolving bias model for the 
galaxy populations. The ratios of the figures of merit (inverse 
areas) are given for various parameter spaces listed in the first 
column. The second column shows the ratios for the Stage IV ex- 
periments alone; the third column includes Stage III information 
for each of them. 

Regarding the number densities used for the galaxy 
populations, these come from selection functions of the sur- 
vey with respect to the intrinsic populations within the de- 
tection limits. In general, target selection is a complicated 
procedure and these numbers represent a sculpted target 
sample not a flux- or volume-limited distribution. We con- 
sider one simple variation in the BigBOSS ELG distribu- 
tion, motivated by the previous subsection where the red- 
shift range was shifted from z = 1 — 2 to z = 0.7 — 1.7. 
Such a shift was found to slightly improve the cosmology 
constraints, and it also reduces the amount of time needed 
to observe the galaxies. If we take advantage of this by now 
looking at a survey plan with four times the number density 
of ELG in the range z = 0.7— 1 (and z = 1 — 1.7 unchanged), 
we find further improvements in determination of 7, wo, w a 
by 2%, 4%, 3% relative to the uniform number density in 
z = 0.7 — 1.7 case of the previous subsection. 

These calculations show that the basic point of ground 
and space surveys being capable of delivering comparable 
cosmology constraints is not very sensitive to these varia- 
tions in the survey design. Detailed experiment design and 
optimization, however, is beyond the scope of this paper. 
We have not consid ered other experimen ts such as the Eu- 
clid space mission (jCimatti et al.l l2009l i. since it includes 
other cosmological probes on a par with the power spec- 
trum me asurement, and 21 cm mapping surveys such as 
SKA (see iPeterson et~ai1 (|2009h : lMoraies fc Wvithd (|2009t )L 
since neutral hydrogen gas measurement techniques and pre- 
cision constraints are not as fully developed. 



5 CONCLUSIONS 

The three-dimensional distribution of large scale structure 
contains information on both the cosmological parameters 
and testing gravity. We have studied the capabilities of next- 
generation power spectrum experiments from the ground, 
BigBOSS, and from space, JDEM-PS, to use the baryon 
acoustic oscillations, power spectrum shape, and redshift 
space distortions to test standard cosmology. 

The main conclusion is that the two experiments could 
achieve comparable constraints. We emphasized the impor- 
tance of including simultaneously the parameters that affect 
growth - the gravitational growth index characterizing de- 
viations from general relativity, the dark energy equation 
of state value and its time variation, and neutrino mass. 
Including these and other cosmological parameters we esti- 
mate the uncertainty on determination of the gravitational 
growth index to be 0.043 for BigBOSS, 0.054 for JDEM-PS, 
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